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摘要 : 太阳 高 能 粒子 (SEP) 事件 是 一 类 重要 的 空间 天 气 灾害 性 事件 ， 其 数值 预报 研究 在 空间 天 
气 预 报 研 究 中 占有 很 重要 的 地 位 。SEP 事件 主要 包括 3 种 类 型 : 与 太阳 耀 斑 爆发 相关 联 的 脉冲 型 
事件 ,与 日 冕 物质 抛射 驱动 的 激 波 相关 联 的 缓 变型 事件 ， 以 及 同时 具有 组 变型 和 脉冲 型 事件 特征 
的 混合 型 事件 。 其中, 缓 变型 SEP 事件 持续 时 间 较 长 并 且 高 能 粒子 强度 较 大 , 对 这 类 事件 的 模拟 
是 当前 研究 的 难点 。 目 前 针对 组 变型 SEP 事件 的 模拟 工作 业已 发 展 了 多 个 理论 和 数值 模型 。 每 
个 模型 都 对 SEP 加 速 和 传播 的 复杂 过 程 作 了 基本 的 假设 , 这 些 模型 的 模拟 结果 能 够 部 分 重 现 观 
测 到 的 SEP 事件 特征 。 而 若 要 提高 预报 SEP 事件 的 能 力 ， 则 需要 将 描述 三 维 日 冕 物质 抛射 驱动 
的 激 波 模 型 与 描述 高 能 粒子 在 行星 际 空间 中 的 加 速 和 传输 的 模型 不合 起 来 , 建立 基于 接近 真实 
的 SEP 加 速 和 传播 的 三 维 太 阳 风 背景 模拟 及 以 激 波 参数 为 输入 的 SEP 模型 。 主要 回顾 了 缓 变型 
SEP 事件 中 粒子 的 加 速 和 传输 方面 的 研究 进展 , 以 及 可 用 于 获取 CME 激 波 传播 参数 的 磁 流 体力 
学 太阳 风 模 型 研究 现状 ; 综述 了 缓 变型 SEP 事件 的 激 波 -粒子 模型 (shock-and-particle model); 
最 后 对 未 来 工作 进行 了 讨论 和 展望 。 

X 键 da: 太阳 高 能 粒子 ; BE: MHD 太阳 风 模式 ; 行星 际 输 运 
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1 引言 


太阳 高 能 粒子 (简称 SEP) 事件 是 能 量 粒子 的 通 量 突然 增强 的 事件 ,一 般 可 以 由 人 造 卫 
星 和 地 面 探测 设备 探测 到 。 美国 国家 海洋 和 大 气管 理 局 (NOAA) 将 SEP 事件 定义 为 能 量 大 
T 10 MeV 的 粒子 在 连续 15 min 以 上 的 时 间 内 数目 超过 10 pfu (cm? - s7? - sr 1) 的 爆发 性 
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5E", SEP 事件 能 够 产生 高 通 量 的 高 能 粒子 分 布 , 这 不 仅 会 影响 航天 器 的 正常 运行 , 而 且 
会 对 航天 员 的 生命 安全 造成 严重 危害 。 在 极端 情况 下 , 能 量 较 高 的 SEP 甚至 可 到 达 地 面 , 产 
生地 面 水 平 事件 (简称 GLE), GLE 是 一 种 极端 SEP 事件 ， 和 通常 的 SEP 事件 (粒子 的 能 
量 范围 通常 在 10 ~ 100 MeV 之 间 ) 相 比 ，GLE 事件 中 粒子 的 能 量 甚至 能 达到 1 GeV", 

SEP 事件 可 分 为 3 28". 缓 变型 SEP 事件 (Gradual SEP events), 脉冲 型 SEP 事件 
(Impulsive SEP events) 以 及 混合 型 SEP 事件 (Mixed SEP events)。 绥 变型 SEP 事件 持续 
HARK, 通常 持续 数 天 , 能 量 粒子 强度 较 大 , 质子 含量 相对 丰富 ; 绝 大 多 数 的 缓 变 型 SEP 
事件 都 伴随 着 日 网 物质 抛射 (简称 CME) 所 驱动 的 激 波 而 发 生 。 脉冲 型 SEP 事件 持续 时 间 
比较 短 , 通常 只 有 数 小 时 ,强度 较 弱 ,并且 电子 含量 相对 丰富 , 一 般 认为 与 浴 斑 过 程 相关 。 另 
Sb, Kallenrode" 指出 , 有 些 SEP. 事件 中 粒子 的 成 分 由 脉冲 型 事件 的 特征 演化 为 缓 变型 事件 
的 特征 ; 即 一 次 SEP 事件 中 , 汉王 和 CME 驱动 的 激 波 对 高 能 粒子 的 加 速 过 程 可 能 都 有 贡 
献 ， 导致 所 谓 的 混合 型 事件 。 

对 SEP 事件 的 实地 观测 内 容 包 括 粒子 的 组 分 、 电 荷 态 、 能 谱 、 强 度 分 布 和 各 向 异性 分 布 
5&^ 7. de 1 给 出 缓 变 型 和 脉冲 型 SEP 事件 各 自 的 特点 ， 从 中 可 以 看 出 , 在 脉冲 型 SEP 事件 
+ 3He/*He 丰 度 的 比值 约 为 1， 而 太阳 风 中 的 3He/^He 丰 度 的 比值 虽然 波动 很 大 但 是 很 少 
超过 0.01, 因此 脉冲 型 SEP 事件 有 时 也 称 富 氨 3 电子 事件 。 另 外 , 脉冲 型 SEP 事件 的 Fe/O 
比例 比 正常 日 园 大 气 中 Fe/O 的 比例 大 10 倍 , 这 些 特征 显示 泡 席 起 源 粒子 在 脉冲 型 事件 中 
占 主导 地 位 。 缓 变型 SEP 事件 中 3He/4He 以 及 Fe/O 的 比例 都 与 日 园 大 气 中 的 成 分 比例 接 
E, 关于 两 类 SEP 事件 的 特征 和 区 别 , Reames” F 1999 年 给 出 了 详细 的 论述 。 


表 1 ”组 变型 和 脉冲 型 SEP 事件 的 特点 


缓 变型 脉冲 型 
粒子 质子 丰富 TEE 
?He/^He « 0.01 约 1 
[Fe/O]/[Fe/O]cornal 约 1 约 10 
持续 时 间 天 小 时 
经 度 分 布 约 180? < 30° 
相关 射电 暴 类 型 II IV II, V 


事件 发 生 率 (太阳 活动 极 大 年 ) 约 10 次 /年 ” 约 1000 次 /年 


SEP 事件 会 对 卫星 和 进行 航天 作业 的 宇航 员 构 成 严重 的 威胁 ,如 能 准确 预报 SEP 事件 ， 
人 们 便 可 以 采取 必要 的 防护 措施 , 保障 卫星 和 星 载 设备 的 安全 ; 同时 也 可 以 为 航天 器 的 故障 
分 析 和 航天 员 的 安全 防护 提供 科学 依据 ; 还 可 以 为 采取 应 变 措施 赢得 宝贵 时 间 ， 尽 可 能 地 降 
低 经 济 损失 。 然 而 SEP 事件 产生 和 发 展 的 物理 机 制 非常 复杂 ,目前 还 没有 被 完全 解释 清楚 ， 
预报 SEP 事件 的 能 力也 不 尽 如 人 意 。 因 此 ,正确 理解 和 恰当 的 动力 学 描述 这 些 机 制 , 对 推 
进 和 提高 空间 天 气 预 报 能 力 是 非常 重要 的 。 相 对 于 脉冲 型 SEP 事件 , 缓 变型 SEP 事件 由 于 
其 高 强度 和 长 持续 时 间 等 特征 造成 的 危害 可 能 更 大 , 由 CME 驱动 的 激 波 加 速 引起 的 缓 变型 
SEP 事件 可 以 造成 近 地 空 间 持续 数 天 的 高 达 数 十 MeV/n 甚至 GeV/n 能 量 的 灾害 性 粒子 辐 
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射 扰动 。 所 以 对 缓 变型 SEP 事件 预报 具有 重要 的 空间 天 气 学 意义 ”, 同时 SEP 事件 也 是 每 
一 个 11 a 太阳 活动 周期 内 空间 高 能 粒子 通 量 的 主要 来 源 ” 。 

对 于 脉冲 型 SEP 事件 的 模拟 , 许多 模型 假设 粒子 的 源 项 是 固定 的 ,通常 假定 源 项 分 布 
在 太阳 表面 或 者 接近 于 太阳 表面 , 然后 求解 粒子 在 行星 际 空间 的 传输 过 程 , 得 到 观测 者 位 置 
的 通 量 前 面 、 各 向 异性 分 布 等 。 不 同 于 脉冲 型 SEP 事件 ,组 变型 SEP 事件 的 粒子 源 是 移动 
的 行星 际 激 波 ， 导 致 缓 变型 事件 的 建 模 难 度 更 大 。 缓 变型 SEP 事件 通常 与 快速 的 CME 相 
X", 这 些 快速 的 CME 加 速 行星 际 激 波 , 使 其 速度 非常 高 , ACE (Advanced Composition 
Explorer) 卫星 在 2003 年 10 月 29 日 和 30 日 探测 到 的 实地 太阳 风速 峰值 高 达 1850 km . s^! 
和 1700 km-s~!, GOES (Geostationary Operational Environmental Satellite) 卫星 记录 到 
行星 际 激 波 加 速 产生 的 一 系列 SEP 事件 , 其 中 还 有 3 次 达到 了 GLE KP”. 快速 的 CME 
驱动 的 行星 际 激 波 对 加 速 粒 子 起 着 主导 作用 , 所 以 提前 几 天 或 几 小 时 预测 大 型 组 变型 SEP 
事件 的 技术 难点 主要 包括 以 下 几 点 ””; (1) CME 在 何 时 、 何 地 发 生 ; (2) 详 述 CME 的 特点 ， 
比如 它 能 和 否 驱 动 激 波 等 ; (3) 确定 CME 驱动 的 激 波 加 速 粒子 到 高 能 量 的 效率 , 同时 被 加 速 的 
粒子 怎样 被 注入 到 行星 际 介 质 ; (4) 被 加 速 的 粒子 和 CME 驱动 的 激 波 在 背景 太阳 风 中 是 怎 
样 传播 的 。 
因此 , 缓 变型 SEP 事件 的 模拟 主要 包括 以 下 3 个 方面 : 

1) 行星 际 激 波 位 形 和 演化 的 描述 ; 

2) 粒子 加 速 和 传输 的 模拟 ; 

3) 随 着 激 波 向 外 传播 ， 激 波 加 速 机 制 和 能 量 粒子 注入 的 连续 检验 。 

总 体 来 讲 , SEP 事件 的 模拟 流程 可 用 图 1 表示 。 对 于 脉冲 型 SEP 事件 , 只 要 将 太阳 高 能 
粒子 的 注入 率 输入 到 粒子 模型 里 即 可 进行 有 效 模 拟 ; 而 对 于 缓 变型 SEP 事件 的 模拟 ， 则 需 
要 首先 运行 磁 流 体力 学 (简称 MHD) 模型 ， 获 取 激 波 面 沿 磁力 线 的 cobpoint" (Connecting 
with the OBserver POINT: 连结 太阳 和 观测 者 的 磁力 线 通 过 激 波 面 上 的 位 置 ) 的 激 波 参 数 
(比如 速度 、 磁 场 )， 作 为 粒子 模型 的 输入 参数 ,然后 通过 求解 传输 方程 得 到 SEP 的 能 谱 、 强 
度 和 各 向 异性 分 布 。 从 图 中 可 以 看 出 , 无 论 是 模拟 缓 变型 还 是 脉冲 型 SEP 事件 ,给 定 比较 准 
确 的 粒子 注入 率 很 关键 。 模拟 脉冲 型 事件 时 ,一般 假设 粒子 源 被 固定 在 太阳 附近 ; 而 缓 变 型 
事件 粒子 源 是 随 太 阳 风 一 起 运动 的 行星 际 激 波 , 高 能 粒子 在 激 波 附近 被 加 速 时 ， 可 以 逃离 激 
波 进 入 太阳 风 并 沿 着 行星 际 磁场 传播 到 地 球 附近 , 而 且 传 输 的 空间 和 时 间 尺 度 比 较 大 , 意 
着 传输 条 件 会 有 显著 改变 , 因此 粒子 注入 率 与 CME 所 驱动 的 激 波 传播 参数 的 准确 获取 有 着 
密 不 可 分 的 关系 。 若 要 建立 更 加 真实 精确 的 SEP 物理 模式 , 则 需要 将 描述 三 维 日 冕 物质 抛射 
驱动 的 激 波 MHD 模型 与 描述 高 能 粒子 在 行星 际 空间 中 的 加 速 和 传输 的 粒子 模型 耦合 起 来 。 
目前 缓 变型 SEP 事件 模拟 的 相关 模型 中 , 一 般 需 要 对 激 波 的 传播 和 SEP 加 速 与 输 运 
的 复杂 过 程 作出 一 些 基本 的 简化 和 假设 “””， 当 模拟 激 波 的 传播 时 ,假设 激 波 面 为 平面 或 球 
面 ; 把 激 波 对 粒子 的 加 速 过 程 当 做 黑匣子 处 理 ， 人 为 给 定 一 个 激 波 加 速 源 项 ; 行星 际 磁场 假 
设 成 Parker 场 等 。 这 些 模型 能 够 得 到 一 些 与 实际 观测 相 吻 合 的 结果 , 但 是 并 不 能 全 面 准确 
地 重 现 SEP 事件 的 观测 特征 , 预报 的 准确 度 较 低 。 因 此 ,融合 接近 真实 的 三 维 MHD 太阳 风 


太阳 风 和 行星 
sewer 


激 波 加 速 
粒子 的 注入 
观测 者 位 置 的 通 量 章 面 


图 1 ” 激 波 -粒子 模型 的 基本 模块 


背景 , 考虑 CME 驱动 的 激 波 的 三 维 结构 ， 以 及 包含 更 真实 丰富 的 物理 内 容 的 激 波 加 速 模式 ， 
是 组 变型 SEP 事件 物理 模式 研究 的 方向 。 

本 文 第 2 章 介绍 SEP 事件 中 粒子 的 加 速 和 传输 方面 的 研究 进展 , 第 3 章 回 顾 可 用 于 获 
取 CME 激 波 传播 参数 的 MHD 太阳 风 模 型 , 第 4 章 综述 缓 变型 SEP 事件 的 激 波 -粒子 模型 ， 
第 5 章 对 未 来 工作 进行 讨论 和 展望 。 


2 SEP 事件 粒子 的 加 速 和 传输 研究 进展 


2.1 Bee SEP 事件 中 粒子 的 加 速 模型 

缓 变型 SEP 事件 的 模拟 需要 考虑 激 波 加 速 的 影响 ， 相 关 模 型 可 以 粗略 地 分 为 两 类 ， 
(1) 不 考虑 具体 的 激 波 加 速 机 制 ”“”， 即 忽略 激 波 对 粒子 的 具体 加 速 过 程 , 在 传输 方程 中 加 
入 一 个 激 波 加 速 源 项 来 描述 粒子 被 运动 激 波 加 速 后 在 激 波 面 附近 的 分 布 函 数 ; 然后 求解 传 
输 方 程 , 得 到 观测 者 处 在 不 同位 置 时 粒子 的 通 量 。 通过 将 模拟 结果 与 观测 结果 作对 比 ， 进 而 
研究 高 能 粒子 注入 效率 函数 随 粒 子 能 量 、 时 间 和 空间 的 变化 ， 以 及 粒子 在 行星 际 的 传播 过 
程 。(2) 发 展 激 波 加 速 理论 , 主要 包括 激 波 漂移 加 速 、 扩 散 激 波 加 速 ”“、 随 机 加 速 ” 等 。 
这 几 种 激 波 加 速 机 制 的 相对 贡献 依赖 于 激 波 的 性 质 ， 激 波 速度 、 压 缩 比 、Alfvén 马赫 数 、 上 
游 磁 场 与 激 波 法 向 的 夹 角 (Opn) 以 及 激 波 附近 存在 的 MHD 扰动 都 能 改变 粒子 加 速效 率 的 相 
对 贡献 大 小 ”。 下 面 分 别 对 这 两 类 工作 进行 介绍 。 
2.1.1 不 考虑 具体 加 速 机 制 的 模型 

Heras 等 人 ”对 缓 变型 SEP 事件 的 空间 演化 进行 了 建 模 研究 , 在 Heras 的 模型 中 首 
次 采用 聚焦 传输 方程 ,用 其 蔡 代 Parker 传输 方程 模拟 SEP 在 行星 际 空间 的 传播 过 程 。 这 个 
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模型 的 主要 特点 如 下 : (1) 在 聚焦 传输 方程 中 加 入 一 个 激 波 加 速 源 项 @, 表征 粒子 被 激 波 加 
速 在 cobpoint 位 置 的 粒子 注入 率 。 运 用 此 传输 方程 ,可 以 很 好 地 再 现 SEP 事件 中 大 而 持久 
的 各 向 异性 分 布 “” 。(2) 对 cobpoint 概念 的 详细 阐述 与 量化 实现 。cobpoint 随 着 激 波 前 沿 改 
变 且 沿 着 连结 飞船 与 太阳 的 磁力 线 运 动 , 通过 追踪 cobpoint 的 运动 轨迹 、 激 波 法 向 、 激 波 
上 下 游 速度 的 变化 规律 ， 可 以 为 不 同位 置 的 飞船 提供 激 波 加 速 粒 子 的 注入 强度 -时 间 变 化 曲 
线 , 这 在 SEP 事件 研究 中 具有 重要 作用 。 模型 假设 激 波 上 游 的 行星 际 磁场 采用 Parker 螺旋 
场 模型 , 粒子 的 注入 发 生 在 cobpoint, 采用 MHD 模型 描述 激 波 从 0.08 ~1 AU (1 AU BK 
阳 到 地 球 的 平均 距离 , 在 2012 年 8 月 于 中 国 北京 举行 的 国际 天 文学 大 会 上 , 该 值 被 定义 为 
149 597870700 m) 的 传播 过 程 ， 从 而 得 到 激 波 参数 在 cobpoint 处 的 演化 规律 。 此 外 , COS 
等 人 中 通过 模拟 cobpoint 沿 着 连结 ACE 与 太阳 磁力 线 的 运动 来 追踪 cobpoint 的 演化 , 并 
且 给 出 了 三 维 运动 图 像 ( 见 图 2)。 由 于 粒子 的 注入 发 生 在 cobpoint ME, 未 来 可 以 提高 获取 
cobpoint 位 置 的 算法 ， 以 期 在 粒子 模型 中 能 够 获取 更 加 准确 的 粒子 注入 率 , 使 模拟 结果 更 加 
接近 真实 。 

在 Heras 等 人 假设 的 基础 上 ，Kallenrode 5& AU “把 激 波 模拟 为 半球 形 的 , 使 其 以 恒定 
速度 沿 径 向 从 太阳 表面 向 外 传播 ; 而 且 激 波 注入 效率 函数 随 径 向 距离 和 方位 角 变 化 , 用 以 体 
现 激 波 加 速效 率 的 时 空 变 化 。 他 们 仅仅 改变 举 焦 尺度 来 修正 粒子 在 激 波 下 游 的 传播 , 但 是 他 
们 没有 考虑 激 波 传 播 经 过 被 扰动 后 的 行星 际 磁场 所 带 来 的 影响 , 这 将 会 导致 不 同 的 结果 ; 他 
们 还 考虑 到 粒子 横越 激 波 的 传播 , 却 没有 考虑 当 粒 子 被 反射 后 或 者 传播 过 程 中 能 量 的 改变 。 

Ng 等 人 "给 出 一 个 高 能 离子 和 Alfvén 波 相互 耦 合演 化 的 模型 ,他们 结合 了 磁 聚 焦 效 
应 、 投 掷 角 扩 散 效 应 、 波 的 放大 效应 和 一 个 运动 的 离子 源 。Ng 等 人 采用 有 限 差分 方法 , 求解 
耦合 的 聚焦 传输 方程 和 描述 高 能 离子 激发 的 等 离子 体 波动 的 强度 演化 方程 。 该 模型 假设 激 波 
以 恒定 速度 沿 径 向 传播 , 参数 化 粒子 注入 效率 随时 间 、 粒 子 刚度 和 径 向 距离 变化 。 该 模型 最 
大 的 优势 是 考虑 了 粒子 在 传输 过 程 中 激发 的 等 离子 体 波 动 对 粒子 散射 的 影响 。 但 是 行星 际 磁 
场 沿 径 向 这 一 假设 过 于 简单 ,不 能 解释 SEP 事件 的 经 度 差异 ; 并 且 激 波 以 恒定 速度 传播 ,是 
缺乏 物理 依据 的 追踪 cobpoint 的 方式 ” 。 

Lario 5& A "^ “进一步 优化 了 Heras 的 模型 , 在 传输 方程 中 加 入 太阳 风 对 流 和 绝热 减速 
效应 。 通过 模拟 不 同 的 SEP. 事件 ,给 出 了 一 个 在 cobpoint 位 置 激 波 加 速 质子 的 注入 率 @ 与 
径 向 速度 比值 VR 的 经 验 关系 式 : lg 8 x VR (VR 来自 MHD 模型 的 输出 , Q 作为 粒子 传输 
模型 的 参数 输入 )。 模 型 成 功 地 再 现 了 一 些 SEP 事件 中 低能 段 (小 于 20 MeV) 质子 通 量 和 各 
ae. Lario 等 人 ”建立 的 MHD 变量 与 粒子 注入 率 的 Q(V R) 关系 , 虽然 只 是 
一 个 半 经 验 关 系 , 却 为 我 们 提供 了 激 波 -粒子 混合 模拟 的 范例 ， 让 我 们 思考 怎样 把 粒子 注入 
率 与 CME 驱动 的 激 波 的 传播 参数 联系 起 来 。 但 是 ,Lario 等 人 提出 的 Q(V R) 关系 式 仅仅 将 
注入 率 与 cobpoint 激 波 的 速度 联系 在 一 起 ; 未 来 可 以 寻找 更 加 真实 的 关系 , 将 MHD 模型 得 
到 的 参量 与 粒子 注入 率 联系 起 来 。 

在 以 上 研究 的 基础 上 , Aran 5 AP ”建立 了 预报 SEP 的 SOLPENCO (SOLar Particle 
ENgineering COde) 模型 。 图 3 给 出 了 SOLPENCO 在 欧洲 空间 天 气门 户 网 站 的 初始 界面 ， 
这 个 操作 工具 包括 两 个 界面 : 需要 输入 不 同 变量 值 (日 心 距离 、 初 始 脉冲 速度 、 观 测 者 角度 位 
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图 2 ”不 同时 间 段 的 三 维 运动 图 像 ”” 


黑色 为 行星 际 磁场 线 , 实心 黄色 圆圈 为 太阳 的 位 置 ,实心 紫色 圆圈 为 ACE, 实心 红色 圆圈 为 cobpoint， 色 标 代 表 相 对 密度 
(p 一 po)/po 的 等 值 面 。(a) t = 30 min; (b) t = 10 h; (c) t = 20 h; (b) t = 30 hi (b) t = 40 h. 
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署 、 平 均 自 由 程 、 是 否 存在 高 度 滑动 、 粒 子 能 量 ) 的 初始 界面 ; 根据 输入 的 参量 , 产生 某 一 给 
定 事件 的 通 量 图 像 或 者 一 组 数值 矩阵 。 此 模型 的 空间 天 气 应 用 详 见 参考 文献 [12]。 
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图 3 SOLPENCO 在 欧洲 空间 天 气门 户 网 站 的 初始 界面 


上 述 模 型 都 没有 涉及 粒子 被 运动 激 波 加 速 的 具体 物理 机 制 , 也 不 考虑 激 波 面 的 MHD 条 
JE. 如 激 波 位 形 、 激 波 强 度 和 速度 演化 如 何 影响 激 波 加 速效 率 等 问题 。 要 对 粒子 加 速 和 注入 
的 复杂 现象 和 机 制 进行 描述 ,只 做 运动 激 波 的 假设 显然 不 充分 。Lario 等 人 "” 提 出 将 激 波 加 
速 粒 子 的 注入 率 的 演化 与 激 波动 力学 性 质 联系 起 来 的 Q(VR) 关系 ,这 个 关系 能 够 量化 注入 
率 等 参数 , 但 是 无 法 解释 激 波 加 速 粒 子 的 物理 机 制 。 
2.1.2 考虑 具体 加 速 机 制 的 模型 

也 有 一 些 学 者 尝试 研究 较为 真实 的 CME 驱动 激 波 附近 粒子 注入 加 速 过 程 。 关 于 激 波 加 
速 的 模拟 工作 主要 集中 在 扩散 激 波 加 速 , 扩散 激 波 加 速 是 一 种 一 阶 费 米 加 速 ” 过 程 。 当 粒子 
多 次 来 回 穿越 激 波 面 后 , 会 被 加 速 到 很 高 的 能 量 , 这 种 加 速 过 程 如 图 4 所 示 。 

下 面 简单 介绍 扩散 激 波 加 速 的 分 析 计 算 , 主要 内 容 见 参考 文献 [132]。 假 设 散 射 在 激 波 
上 下 游 都 起 作用 ， 引 起 粒子 沿 着 磁力 线 以 及 垂直 磁力 线 的 方向 扩散 ; 等 离子 体 分 布 范围 从 
x= —oo 到 zz 二 十 co, 无 限 湾 的 平面 激 波 处 于 z = 0 处 , 粒子 的 速度 远大 于 背景 流体 速度 ; 
那么 描述 粒子 输 运 的 稳 态 玻 尔 兹 曼 方 程 可 以 写成 扩散 对 流 方程 的 形式 ” : 
ED zu IOF Dp fm»). (1 


x 


2 [urten - 


4 ”扩散 激 波 加 速 过 程 示意 图 
垂直 阴影 线 为 激 波 面 , 圆 点 表示 散射 元 ， 带 有 稍 头 的 实 线 表 示 加 速 粒子 的 路 径 。 


其 中 , /是 粒子 的 分 布 函数 , k 是 垂直 于 激 波 面 的 扩散 系数 , wu 是 流体 速度 , p 是 粒子 的 动量 。 
由 式 (1) 可 以 推出 上 游 和 下 游 无 穷 远 边界 处 的 分 布 函数 fi 和 fo 的 关系 表达 式 , 为: 


fip) -tw [ari mna b=, (2) 


uj 一 U2 


其 中 , u 和 wz 分 别 代表 上 游 和 下 游 流 体 速度 ， 户 是 激 波 下 游 无 穷 远 处 的 分 布 函数 ， 广 是 激 
波 上 游 无 穷 远 处 的 分 布 函数 。 这 里 不 给 出 此 公式 的 推导 过 程 , 详细 的 推导 见 参考 文献 [19]。 

假设 注入 粒子 能 谱 是 单一 能 量 的 , 即 fa (p) ox 5(p — po)， 那 么 从 式 (2) 可 以 得 出 , 当 
p > po Hl. hp x p EREN. RRM b BEAL y = /ua RE, 根据 马赫 数 
M 二 w/vai (var 为 激 波 上 游 Alfvén 速度 ), 可 以 得 到 : 

3r 3(P- 1M? (T +1)M? 
^Cr-17 £*«Mi-1)" rTM (3) 

其 中 , P= 5/3 为 气体 的 绝热 指数 。 当 M? — oo, 强 激 波 条 件 下 压缩 比 r — 4, BEIM 
b — 4， 而 观测 到 的 字 宙 线 能 谱 为 5 = 4.7, 因此 这 种 加 速 机 制 常 被 用 于 解释 字 宙 线 的 起 源 

扩散 激 波 加 速 理论 基于 Parker PREM: 


2 


十: apV ， Uo + oss (Duo E ， (4 
Rp, 了 是 粒子 分 布 函 数 ， U 为 背景 流 场 的 速度 ,KK 和 Do 分 别 为 空间 扩散 系数 张 量 和 动量 
PRAM, Va 是 漂移 速度 。Parker 传输 方程 能 融合 运动 电场 导致 的 漂移 加 速 ， 粒 子 穿越 激 
波 在 上 下 游 介质 连续 散射 的 一 阶 费 米 加 速 ， 由 等 离子 体 湛 动 引起 的 随机 激 波 加 速 等 机 制 ， 可 
以 成 功 解释 CME 驱动 的 激 波 加 速 产生 的 单 贤 律 谱 能 谱 特 征 以 及 其 他 观测 现象 , 但 前 提 必 须 
假设 粒子 投掷 角 分 布 各 向 同性 。 

对 于 扩散 激 波 加 速 的 模拟 工作 主要 有 以 下 几 点 所 :; (1) 对 于 CME 驱动 的 激 波 ， 由 于 
其 加 速 粒子 的 时 间 和 能 量 有 限 ， 因 此 粒子 要 获得 加 速 ， 需 要 在 激 波 附 近 停 留 足够 长 的 时 
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EPST, (2) 逃逸 的 高 能 粒子 在 激 波 上 游 可 以 激发 并 放大 Alfvén W, 这 些 被 放大 的 Alfvén 
波 能 够 有 效 地 把 粒子 保持 在 激 波 附近 。Lee” 提出 了 一 个 与 磁 流 体力 学 波 相 关 的 扩散 激 波 
加 速 粒 子 的 自 洽 理论 , 研究 表明 Alfvón 波 的 强度 和 加 速 粒子 分 布 函数 的 各 向 异性 程度 相关 ， 
由 此 可 以 得 出 Alfvén 波 的 放大 过 程 和 粒子 的 加 速 过 程 存在 耦合 。 此 外 , Ng 等 人 "在 模拟 组 
变型 SEP 事件 时 建立 了 激 波 加 速 高 能 粒子 与 行星 际 Alfvén 波 相 耦合 的 理想 模型 ,模拟 结果 
表明 粒子 激发 的 Alfvén 波 对 它们 的 传输 过 程 具有 重要 作用 。(3) 较为 真实 的 CME 驱动 的 激 
波 并 不 是 简单 的 平面 或 球面 激 波 , 而 是 具有 有 限 经 度 宽 的 复杂 三 维 结构 ,其 演化 也 受 三 维 太 
阳 风 结构 的 影响 。Sokolov 等 人 "通过 克 合 一 个 三 维 的 MAD 模型 和 粒子 加 速 的 场 线 平流 模 
型 (field line advection model for particle acceleration), 得 到 一 个 由 较为 真实 的 CME 驱动 
的 激 波 加 速 粒子 的 模型 。 该 模型 能 得 到 一 个 较为 真实 的 CME 和 激 波 结构 ， 并 模拟 粒子 在 距 
离 太 阳 4 ~ 30 R 范围 内 的 加 速 过 程 , 结果 发 现 扩散 加 速 能 够 将 粒子 加 速 到 几 百 MeV。 其 模 
拟 结果 与 1998 年 5 月 2 日 发 生 的 SEP 事件 的 一 些 观测 结果 吻合 较 好 。 

Zank ZAE TA “HAR” A (Onion Shell Model), 发 展 了 基于 扩散 激 波 
IR MAE St. “RE” 模型 是 粒子 加 速 和 波 激发 过 程 的 耦合 , 可 近似 为 一 个 激 
波 掠 过 的 同心 球 层 , 而 忽略 星际 介质 可 能 带 来 的 不 对 称 性 。 该 模型 假设 行星 际 磁场 为 Parker 
场 , 爆炸 激 波 在 球 对 称 太阳 风 中 传播 , 采用 一 维 的 ZEUS MHD 模式 模拟 激 波 产生 后 在 内 日 
球 层 中 的 传播 过 程 ， 激 波 参 量 如 强度 、 马 赫 数 等 和 太阳 风 的 演化 由 流体 力学 数值 模拟 确定 。 
在 准 平行 激 波 条 件 下 , 粒子 的 加 速 和 Alfvén 波 的 非 线 性 增长 耦合 在 一 起 ,粒子 在 传播 、 演 化 
的 行星 际 激 波 处 得 以 加 速 。 模 型 还 包括 激 波 加 速 后 的 粒子 逃逸 到 上 游 激发 Alfvén 波 的 过 程 。 
模拟 结果 显示 ,行星 际 激 波 加 速 粒子 能 达到 的 最 大 能 量 由 逐渐 降低 的 激 波 速度 、 逐 渐 减弱 的 
行星 际 磁场 和 激 波 传播 时 间 共 同 决定 。 

Rice 等 人 中 进一步 优化 了 Zank 等 人 的 模型 ， 利 用 波 增长 方程 的 稳 态 解 自治 地 计算 并 得 
到 加 速 粒子 激发 的 上 游 波 强度 , 这 样 空间 扩散 系数 也 可 以 通过 自治 地 计算 获得 , 空间 扩散 系 
数 与 动量 相关 并 且 决 定 粒子 的 加 速 和 随后 的 演化 。 该 模型 适用 于 任意 强度 的 激 波 ,明显 地 改 
进 了 Zank 等 人 通常 只 适用 于 强 激 波 的 模型 。 

李刚 等 人 "对 Zank 和 Rice 等 人 建立 的 模型 做 了 较 大 改进 , 采用 蒙特 卡 洛 方法 模拟 
粒子 经 激 波 加 速 后 ,在 距 激 波 较 远 处 逃逸 后 在 Paker 螺旋 行星 际 磁场 的 传播 过 程 , 该 模型 可 
以 求 得 粒子 往返 穿越 1 AU 的 次 数 , 继而 可 以 较为 准确 地 计算 出 事件 中 粒子 的 强度 分 布 和 能 
dE. 此外， 李刚 等 人 还 运用 该 方法 研究 了 行星 际 满 流 对 粒子 输 运 的 影响 , 模型 给 出 了 不 同 平 
均 自由 程 、 不 同 能 量 粒 子 的 强度 曲线 随时 间 的 变化 , 并 且 理 论 模 型 可 以 很 好 地 解释 ACE 和 
WIND 飞船 的 观测 结果 。 在 此 模型 的 基础 上 ,李刚 等 人 “计算 了 同时 有 炮 斑 和 CME 驱动 激 
波 加 速 时 的 时 间 强 度 剖面 ,结果 显示 : 时 间 强 度 剖面 经 历 一 段 快 速 增长 ， 即 呈现 出 只 有 耀 斑 
加 速 时 的 特征 ， 然 后 趋 于 平稳 状态 ， 即 呈现 出 只 有 激 波 加 速 时 的 特征 。 这 与 Cane 等 人 "” 提 
出 的 一 些 大 SEP 事件 中 某 些 高 能 粒子 的 时 间 强 度 曲线 会 有 两 个 峰值 的 观测 结果 对 应 。 同 时 ， 
模型 还 预测 了 处 于 不 同 经 度 的 观测 者 所 观测 到 的 时 间 强 度 曲 线 ( 见 图 5)。 

后 来 该 模型 发 展 成 为 一 套 较 成 熟 的 日 球 层 粒子 传输 和 加 速 模型 ， 称 为 PATH (Particle 
Acceleration and Transport in the Heliosphere) guu"^ PATH 模型 采用 的 是 基于 Parker 
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图 5 ”模型 预测 的 处 于 不 同 经 度 处 的 观测 者 观测 到 的 时 间 强 度 曲线 
不 同 颜色 的 线 型 代表 不 同 能 量 的 粒子 的 时 间 强 度 。(b) 表示 观测 者 只 与 耀 斑 相连 ; (c) 表示 观测 者 只 与 CME 驱动 的 激 波 相 
连 ; (d) dero BN E SRAM CME 驱动 的 激 波 都 相连 。 


传输 理论 (基于 Parker 传输 方程 ) 的 扩散 激 波 加 速 机 制 , 该 模型 成 功 地 模拟 了 几 例 大 的 缓 变 
型 SEP 事件 的 粒子 相对 强度 的 时 间 剖 面 和 能 谱 ， 计 算 结 果 和 观测 数据 吻合 较 好 。PATH 模 
型 的 成 功 之 处 是 融合 了 准 平行 激 波 的 一 阶 费 米 加 速 机 制 、 激 波 传播 过 程 中 强度 速度 演化 、 太 
阳 风 满 流 对 激 波 加 速效 率 的 影响 , 包含 更 真实 丰富 的 物理 内 容 。 

李刚 等 人 "在 PATH 模型 中 加 入 任意 的 磁场 上 游 与 激 波 法 向 的 夹 角 bo， 用 以 研究 斜 激 
波 上 的 粒子 加 速 , 并 且 考 虑 go。 对 注入 能 量 和 注入 效率 ,平行 扩散 系数 I HHA RAM 
ki ， 以 及 总 的 扩散 系数 的 影响 。 李刚 等 人 自治 地 同时 求解 粒子 和 波 的 输 运 方程 得 到 平行 扩 
WAM kp 计算 垂直 扩散 系数 时 采用 非 线性 引导 理论 。 模拟 结果 显示 , 粒子 注入 效率 随 
Bon 的 增加 迅速 下 降 , 接近 于 太阳 时 (r < 0.2 AU); 准 平行 激 波 比 准 垂直 激 波 能 加 速 粒子 到 
更 高 的 能 量 , 原因 在 于 质子 流 激发 的 Alfvén 波 在 粒子 加 速 中 发 挥 主要 作用 。 准 平行 激 波 加 速 
的 粒子 能 得 到 的 最 大 能 量 随 着 日 心 距离 r 比 准 垂直 激 波 下 降 得 更 迅速 , 而 且 准 平行 激 波 上 游 
波 强度 随 着 距离 下 降 得 更 快 ; 因此 , 随 着 日 心 距 离 的 增加 , 准 垂直 激 波 能 将 粒子 加 速 到 更 高 
的 能 量 。 
2.2 ”各 向 异性 激 波 加 速 粒子 模型 

ATHINE, 早期 人 为 给 定 激 波 加 速 源 项 的 处 理 方式 显然 无 法 精确 描述 激 波 加 速 的 高 
能 粒子 源 分 布 , 而 现 有 的 高 能 粒子 在 CME 驱动 激 波 的 加 速 模型 采用 的 都 是 基于 Parker 传 
输 理 论 的 扩散 激 波 加 速 (比如 PATH 模型 )。Parker 传输 方程 可 以 成 功 解释 某 些 缓 变型 或 者 
混合 型 SEP 事件 的 观测 特征 尤其 是 能 谱 特 征 。 但 是 Parker 传输 理论 是 基于 扩散 框架 下 的 ， 
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以 一 个 基本 假设 为 前 提 ， 即 粒子 投掷 角 分 布 是 各 向 同性 的 ; 显然 ，Parker 传输 理论 不 能 应 用 
于 各 向 异性 粒子 环境 下 的 激 波 加 速 问题 。 而 且 最 近 的 一 些 SEP 事件 的 观测 结果 与 扩散 激 波 
加 速 理论 不 完全 吻合 。 如 Terasawa SA” 发 现 加 速效 率 与 Alfvén 马赫 数 相关 , 而 Alfvén 马 
赫 数 取决 于 上 游 磁场 强度 ; 在 扩散 激 波 加 速 理 论 下 , 特别 是 平行 激 波 时 ， 加 速效 率 只 取决 于 
激 波 压缩 比 , 而 与 上 游 参 数 无 关 。 

为 了 研究 磁 聚 焦 效 应 、 由 汕 动 磁场 导致 的 投掷 角 扩 散 传输 效应 以 及 绝热 冷却 效应 对 SEP 
传输 机 制 的 影响 , 秦刚 等 人 "” 运用 复合 磁 灌流 模型 ,用 数值 模拟 的 方法 计算 了 带电 粒子 在 
Parker 场 的 运动 轨迹 , 并 将 模拟 结果 与 解析 公式 的 计算 结果 进行 比较 。 相对 于 各 向 同性 理 
it, 模拟 结果 与 各 向 异性 绝热 冷却 效应 理论 更 加 吻合 , 并 且 , 模拟 结果 与 磁场 聚焦 和 投掷 角 
散射 理论 得 出 的 结果 基本 一 致 。 

近年 来 , 一 些 研究 者 ”开始 用 聚焦 传输 方程 来 描述 粒子 在 激 波 附近 的 传输 和 加 速 过 
#2, 聚焦 传输 方程 的 本 质 是 Fokker-Plank 方程 ,用 于 描述 粒子 分 布 函数 S 在 相 空 间 的 演化 ， 
其 形式 如 下 : 


Of Jo vf- l Doy OL 
ót =V -kı -Vf —(vub-- Vaw): VF + au Out 
(ene pH) y app. of 


|; Vaw — bb: V Vew) + bb: vv. pof 4 
2 Op 
Ep, p 为 投掷 角 的 余弦 值 , k, JÉSEELU MAR, Du, ARMA BAM, Lp = (b-Vln B)? 
为 磁场 察 焦 尺度 , b 是 沿 着 磁场 方向 的 单位 矢量 , v 是 粒子 的 速度 ,Vs 是 太阳 风 的 速度 。 聚 焦 
传输 方程 唯一 的 假设 是 基于 小 扰动 的 准 线性 近似 。 相 对 于 Parker 传输 方程 ,该 方程 不 仅 融合 
了 费 米 加 速 、 激 波 漂移 加 速 、 横 越 激 波 电场 引起 的 加 速 等 , 而 且 保 留 了 投掷 角 信息 , 不 用 假设 
粒子 各 向 同性 ,所 以 基于 聚焦 传输 方程 的 激 波 加 速 模式 一 般 称 为 各 向 异性 激 波 加 速 模 式 呈 。 

左 平 兵 等 人 "采用 前 向 随机 微分 方法 的 各 向 异性 激 波 加 速 模 式 , 在 连续 单一 能 量 粒子 注 
入 条 件 下 研究 不 同类 型 激 波 (平行 激 波 、 斜 激 波及 准 垂直 激 波 ) 的 加 速 能 谱 ， 比 较 了 各 向 异性 
加 速 和 基于 Parker 传输 方程 的 各 向 同性 加 速 的 异同 点 。 对 于 能 量 较 高 的 粒子 ， 能 谱 接 近 于 
ARE, 这 与 标准 的 扩散 激 波 加 速 理论 一 致 ， 对 于 和 斜 激 波 或 者 准 垂直 激 波 , 能 谱 的 高 能 成 分 
表现 出 双 窜 律 分布 。 当 粒子 注入 源 项 相同 时 ， 由 聚焦 传输 方程 得 出 的 能 谱 强 度 不 同 于 标准 的 
扩散 激 波 加 速 理论 得 出 的 能 谱 强度 。 图 6 (a) 给 出 了 基于 各 向 异性 激 波 加 速 模式 下 平行 激 波 
加 速 高 能 粒子 下 游 的 能 谱 ， 实 线 代 表 运 用 前 向 随机 微分 方法 得 出 的 模拟 结果 ,点 划 线 代表 一 
HE Paker 传输 方程 的 解析 解 。 从 图 6 (a) 中 可 以 看 出 ， 较 高 能 段 模拟 得 到 的 寡 律 分 布 的 强度 
小 于 标准 的 扩散 激 波 加 速 理论 的 值 ， 这 意味 着 很 少 的 粒子 能 被 加 速 到 较 高 的 能 量 。 此 外 , 左 
平 兵 等 人 "还 用 此 模式 研究 了 激 波 加 速 的 高 能 粒子 的 空间 分 布 和 各 向 异性 , 模拟 发 现 激 波 上 
游 的 高 能 粒子 处 于 强 各 向 异性 , 但 是 下 游 却 接近 于 各 向 同性 , 并且 激 波 附近 出 现 强度 尖峰 结 
构 。 这 些 模拟 结果 可 以 很 好 地 解释 Voyager 1&2 飞船 在 终止 激 波 附 近 探测 到 的 终止 激 波 粒子 
(能 量 在 几 keV 至 20 MeV 之 间 的 高 能 粒子 ) 的 异常 观测 特征 “, 包括 : (1) 在 终止 激 波 上 游 
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区 域 , 终止 激 波 粒子 具有 很 强 的 投掷 角 各 向 异性 ,并 且 随 时 间 粒 子 各 向 异性 变化 较 快 ， 而 在 
终止 激 波 下 游 ,终止 激 波 粒子 则 显示 出 稳定 的 投掷 角 各 向 同性 分 布 ; (2) 在 终止 激 波 附近 出 
现 强 各 向 异性 强度 尖峰 结构 ; (3) 在 终止 激 波 上 下 游 都 能 观测 到 多 容 律 谱 能 谱 分 布 , 而 不 是 
经 常 观测 到 的 宇宙 线 单 寡 律 谱 分 布 。 
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图 6 (a) 平行 激 波 下 游 的 能 谱 分 布 ; (b) 粒子 动量 与 投掷 角 随 注入 粒子 投掷 角 的 变化 “ 
(a) 实 线 代表 运用 前 向 随机 微分 方法 得 出 的 模拟 结果 ， 点 划 线 代表 一 维 Paker 传输 方程 的 解析 解 。 


2.8 ”垂直 扩散 对 SEP 在 行星 际 空间 传播 的 影响 

在 行星 际 空间 ,存在 着 由 太阳 风 的 波动 和 满 流 引 起 的 随机 小 尺度 不 规则 磁场 ， 由 于 这 些 
随机 扰动 , SEP 在 行星 际 空间 传播 时 运动 轨道 发 生 散 射 ， 导 致 扩散 运动 。 在 前 面 提 到 的 模拟 
SEP 传播 的 模型 中 , 普遍 假定 粒子 只 沿 着 行星 际 磁力 线 方向 扩散 , 而 不 考虑 垂直 于 磁力 线 的 
扩散 过 程 。 在 准 线性 理论 中 , 垂直 扩散 是 由 于 磁力 线 的 随机 行走 引起 , 该 理论 分 别 求 出 平 
行 扩散 系数 和 垂直 扩散 系数 的 表达 式 , 结果 表明 垂直 扩散 系数 相对 于 平行 扩散 系数 可 以 忽略 
Tib. 但 是 Dwyer"“、 张 明 等 人 "通过 分 析 卫 星 的 观测 数据 , 发 现在 一 些 SEP 事件 中 垂直 
扩散 系数 和 平行 扩散 系数 的 比值 可 以 接近 甚至 超过 1。Matthaeus"“ 提 出 非 线性 引导 中 心理 
i£ ( NLGC ) 并 用 其 描述 粒子 在 磁场 中 的 垂直 扩散 , 这 是 首 个 和 数值 模拟 吻合 得 较 好 的 理论 
模型 ， 对 后 来 关于 垂直 扩散 系数 的 研究 具有 重要 的 参考 价值 。 

ak Ag AST SEP 在 三 维 行星 际 空间 的 传播 过 程 ， 并且 对 比 了 不 加 入 垂直 扩散 和 
加 入 垂直 扩散 后 的 模拟 结果 。 他 们 发 现 , 加 入 垂直 扩散 后 , 通 量 分 布 更 加 均匀 ， 训 变相 期 间 
SEP 各 方向 的 通 量 比值 更 加 接近 于 1， 这 是 因为 垂直 扩散 能 很 大 程度 地 抹 掉 行星 际 空间 中 
SEP 通 量 分 布 的 不 均匀 性 。2001 年 9 月 24 日 在 1 AU WE, Wind 飞船 观测 到 在 SEP 事 
件 的 初始 时 刻 , 粒子 还 处 于 各 向 异性 分 布 时 ,朝向 太阳 运动 的 粒子 异常 增加 的 情况 。 秦刚 等 
人 "从 垂直 扩散 的 角度 解释 了 这 一 异常 情况 。 

汪洋 等 人 ”研究 了 垂直 扩散 对 被 CME 驱动 的 行星 际 激 波 加 速 的 高 能 粒子 的 影响 ， 以 及 
高 能 粒子 在 三 维 Parker 场 中 的 传播 过 程 ; 计算 了 当 激 波 传 向 不 同 的 经 度 和 纬度 方向 时 , 位 于 
1 AU 的 观测 者 测量 到 的 粒子 通 量 和 各 向 异性 剖面 随时 间 变 化 的 情况 。 模 型 假设 CME 驱动 
的 行星 际 激 波 是 一 个 运动 的 SEP YR, 源 的 强度 随时 间 和 空间 的 变化 过 程 由 经 验 模型 给 定 , 不 
考虑 粒子 的 加 速 过 程 。 模拟 结果 显示 ,， 当 加 入 垂直 扩散 后 , 通 量 的 初始 时 间 、 峰值 时 间 、 峰值 
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强度 与 事件 的 持续 时 间 都 会 受到 显著 的 影响 ; 并 且 在 考虑 垂直 扩散 效应 后 ,高 能 粒子 能 够 横 
越 磁力 线 传播 。 观测 者 能 够 在 其 磁力 线 与 激 波 面 连结 起 来 之 前 观测 到 高 能 粒子 ， 当 观测 者 所 
在 的 磁力 线 与 激 波 面 断 开 以 后 , 同样 还 有 粒子 不 断 地 注入 或 者 离开 观测 者 所 在 的 磁力 线 。 图 
7 给 出 有 垂直 扩散 和 没有 垂直 扩散 的 情况 下 ，5 MeV 粒子 的 通 量 和 各 向 异性 随时 间 变 化 的 曲 
线 。 对 于 事件 W70 , 有 垂直 扩散 的 情况 下 , 通 量 衰减 得 更 慢 一 些 ; 对 于 N40 事件 , 观测 者 的 
磁力 线 一 直 与 激 波 面 是 断 开 的 , 所 有 的 粒子 只 能 通过 垂直 扩散 才能 到 达观 测 者 。 

秦刚 等 人 "用 与 汪洋 等 人 "相似 的 模型 (模型 不 同 之 处 在 于 投掷 角 扩散 系数 D, ME 
直 扩散 系数 k 的 公式 不 一 样 , 汪洋 等 人 采用 的 是 基于 量 纲 分 析 的 较 简单 的 经 验 公式 , TH 
刚 等 人 采用 的 是 基于 满 流 理论 的 公式 ) 研究 了 卫星 处 于 不 同 经 纬度 和 径 向 距离 时 ， 形 成 著 水 
His (BIT AR SEP 事件 的 衰退 期 ， 分 开 较 远 的 飞船 观测 到 高 能 粒子 的 通 量 几乎 一 致 ， 
通 量 的 差别 在 2 ~ 3 倍 左右 , 并 且 通 量 随 时 间 线 性 演化 的 现象 ) 所 需 的 垂直 扩散 系数 的 大 小 。 
模拟 结果 显示 ,， 当 垂直 扩散 系数 增加 时 ,粒子 通 量 在 空间 分 布 的 不 均匀 性 将 会 变 小 , 荤 水 池 
现象 变 得 更 加 明显 。 然而 , 如 果 激 波 加 速 强度 随 径 向 距离 + 减 小 得 比 r725 慢 时 , 基于 非 线性 
引导 理论 对 垂直 扩散 系数 的 限制 , 那么 车 水 池 现 象 可 能 会 消失 。 因 此 , 对 缓 变型 SEP 事件 的 
蔷 水 池 现 象 的 观测 可 以 被 用 来 定性 地 检测 粒子 扩散 和 激 波 加 速 理 论 。 

从 上 面 汪洋 等 人 的 模拟 工作 中 可 以 得 出 ,即便 垂直 扩散 比 平行 扩散 弱 得 多 ,垂直 扩散 仍 
然 在 粒子 的 传播 过 程 中 起 着 不 可 忽视 的 作用 。 垂 直 扩 散 可 以 用 于 解释 许多 观测 现象 , 特别 是 
当 观 测 者 所 在 的 位 置 没有 通过 行星 际 磁 力 线 与 源 区 直接 相连 时 ， 垂 直 扩散 是 粒子 到 达观 测 点 
的 最 可 能 途径 。 

2.4 SEP 传输 方程 的 求解 

研究 高 能 粒子 在 行星 际 空间 中 的 传输 时 需要 求解 聚焦 传输 方程 。 一 些 学 者 ””" 采用 
有 限 差分 方法 对 方程 进行 离散 求解 。Lario 等 人 ”利用 时 间 分 步 法 将 传输 方程 ( 见 式 (7)) 分 
裂 成 4 个 一 维 方程 (jw 一 ,7 一 ,p 一 , Kz 一 FFE)» p 方程 解释 了 一 半 的 由 产生 的 效应 ; r 一 和 
p- 这 两 个 方程 分 别 解释 ” 和 p 的 变化 ; 最 后 uo 方程 解释 了 由 人 产生 的 另 一 半 的 效应 。 这 4 
个 一 维 方程 在 每 1/4 的 时 间 步 用 有 限 差 分 格式 求解 (根据 每 个 方程 的 特点 选择 具体 的 差分 格 
式 )。 对 于 维 数 较 低 的 传输 方程 , 利用 有 限 差分 方法 可 以 得 到 较 高 精度 的 数值 模拟 结果 ; 但 当 
传输 方程 维 数 较 高 时 , 数值 计算 的 效率 就 会 很 低 。 

近年 来 ,随机 微分 方程 方法 被 一 些 学 者 用 来 求解 聚焦 传输 方程 ， 因 为 聚焦 传输 方程 的 本 
质 是 Fokker-Plank 方程 , 而 在 数学 上 也 可 以 证 明 随 机 微分 方程 和 Fokker-Planck 方程 是 相互 
等 价 的 。 使 用 这 种 方法 , 聚焦 传输 方程 被 分 解 成 一 个 粒子 的 两 项 单 粒 子 运动 : 一 项 是 确定 
项 , 一 项 是 随机 项 。 随机 项 描述 了 粒子 的 一 个 扩散 过 程 ”。 李 刚 等 人 "利用 把 传输 方程 分 
解 成 一 个 高 能 粒子 的 两 项 单 粒子 运动 , 提出 用 蒙特 卡 洛 方法 来 研究 带电 粒子 在 缓 变型 SEP 
事件 中 的 传播 。 模型 中 确定 项 的 部 分 描述 粒子 在 Parker 螺旋 磁场 下 的 运动 ,随机 部 分 描述 了 
粒子 和 太阳 风 满 流 相互 作用 的 扩散 过 程 。 利用 蒙特 卡 洛 方法 模拟 高 能 粒子 在 行星 际 空间 中 传 
输 过 程 的 好 处 是 程序 相对 比较 简单 ,特别 是 在 研究 脉冲 型 SEP 事件 时 ， 如 果 不 需要 考虑 具 
体 的 加 速 过 程 ,可 以 通过 蒙特 卡 洛 方法 模拟 得 到 在 任意 空间 和 时 间 的 能 谱 "。 张 明 等 人 "" 开 
发 了 一 种 计算 SEP 在 三 维 日 球 层 磁场 的 传播 模式 ， 该 模式 采用 后 向 随机 微分 方程 解 聚焦 传 


0 2 4 6 8 100 2 1 6 8 10 
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图 7 ”在 有 垂直 扩散 ( 实 线 ) 与 没有 垂直 扩散 (虚线 ) 的 情况 下 , 对 比 5 MeV 质子 的 通 量 和 各 向 异性 
虚 垂 线 标示 激 波 通过 1 AU 的 时 刻 ; 两 个 短 垂 线 标示 观测 者 所 在 的 磁力 线 与 激 波 面相 连 的 时 间 段 ; 不 同 格子 里 面 的 图 对 应 于 
激 波 被 指向 不 同 的 方向 的 情况 ; 源 的 强度 随 着 径 向 和 人 角 距 衰减 ，S ~ (r/r.) ? exp( 一 |8|/15°); RRA ^ E" 的 曲线 显示 的 是 
cobpoint 效率 随时 间 的 变化 。 


输 方 程 。 模 拟 结果 显示 了 观测 者 处 于 行星 际 空间 不 同 的 经 度 、 纬度 以 及 径 向 距离 时 , SEP 通 
量 和 各 向 异性 随时 间 的 演化 过 程 。 该 模式 的 优势 是 日 球 层 磁场 并 不 局 限于 Parker 螺旋 磁场 ， 
任意 位 形 的 三 维 日 球 层 磁场 都 可 以 作为 模式 参量 , 而 且 容 易 扩展 到 高 维 , 容易 实现 并 行 计算 ， 
数值 误差 小 。 图 8 比较 了 差分 方法 和 随机 微分 方程 方法 的 计算 结果 , 模型 中 忽略 了 绝热 冷却 
和 垂直 扩散 ,最 后 得 到 两 者 的 计算 差别 在 5% UA, 但 是 , 随机 微分 方程 方法 数值 模拟 在 计 
算 时 间 上 有 明显 的 优势 。 


3 MHD 太阳风 模型 


SEP 事件 中 的 高 能 粒子 在 不 同 的 时 间 尺 度 内 被 炮 斑 过 程 或 者 CME 驱动 的 激 波 加 速 , 并 
且 在 被 扰动 后 的 行星 际 磁 场 和 太阳 风 中 传输 ,这些 过 程 都 紧 紧 依赖 于 太阳 风 背 景 场 。 从 图 1 
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[57] 


图 8 ”差分 方法 和 随机 微分 方程 方法 模拟 高 能 粒子 的 全 向 通 量 和 各 向 异性 随时 间 变 化 的 对 比 图 


中 也 可 以 看 出 , 模拟 缓 变型 SEP 事件 时 ，CME 及 其 驱动 激 波 的 传播 参数 与 粒子 注入 率 有 着 
密 不 可 分 的 关系 。 因 此 ， 获 取 更 加 接近 物理 真实 的 背景 场 是 模拟 SEP. 事件 不 可 或 缺 的 部 分 。 
目前 对 背景 场 的 模拟 主要 采用 MHD 模拟 ,下面 简要 介绍 可 用 于 获取 CME 激 波 传播 参数 的 
MHD 太阳 风 模型 的 研究 结果 ,以 期 将 包含 更 真实 丰富 物理 内 容 的 太阳 风 背 景 加 入 到 SEP 物 
理 模式 中 。 

描述 日 冕 /行星 际 激 波 传播 的 模型 中 ， 有 较 简单 将 激 波 模拟 为 半球 形 、 并 以 恒定 速度 沿 
径 向 从 太阳 表面 向 外 传播 的 模型 拉 T, 也 有 复杂 的 MHD 模型 。 基 于 MHD 的 三 维 数值 模拟 
的 模型 中 , 根据 不 同 的 太阳 风速 度 、 密 度 、 温 度 等 边界 条 件 求解 磁 流 体力 学 方程 组 而 获得 稳 
态 的 背景 太阳 风 分 布 。 在 此 基础 上 ， 以 不 同 的 方法 加 入 CME 或 激 波 等 扰动 , 对 CME 和 激 波 
的 传播 过 程 等 进行 分 析 " 。 冯 学 尚 等 人 "从 物理 、 格式 和 计算 方面 对 日 园 行 星际 过 程 的 三 维 
数值 模拟 工作 中 要 注意 的 问题 进行 了 较 完整 的 总 结 , 下 面 对 一 些 模 型 进行 简单 的 介绍 。 
3.1 ”空间 天 气 框架 模型 

空间 天 气 框架 模型 (SWMF) 由 美国 密 西 根 大 学 牵头 的 空间 环境 建 模 中 心 (CSEM) 及 
其 他 10 多 家 合作 单位 共同 开发 。SWMF 包括 日 晃 模 型 、 太 阳 爆 发 事件 模型 、 内 日 球 层 模 
型 以 及 其 他 区 域 模型 。 各 个 模块 可 以 独立 运行 ,也 可 通过 标准 接口 耦合 起 来 "" 。SWMF 的 
日 冕 模型 、 内 日 球 层 模型 及 其 他 区 域 模型 均 采 用 Roe MAR” (Block Adaptive Tree 
Solar-wind Roe-type Upwind Scheme: 块 自 适应 树 型 太阳 风 Roe 迎风 格式 ，BATS-R-US)， 
采用 块 自 适应 树 型 结构 网 格 , 可 以 根据 解 的 物理 性 质 对 计算 网 格 进行 加 密 或 放 粗 。T6th 等 
AS SHAT SWMF 的 整体 结构 及 其 在 各 方面 的 最 新 发 展 。 
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Manchester 5& A" 7? Hj BATS-R-US 代码 建立 一 个 三 维 MHD 数值 模型 来 描述 CME 
AABE 1 AU 的 演化 过 程 ; 在 边界 条 件 的 确定 中 , 此 模型 固定 底部 开放 场 区 的 各 物理 
量 , 然后 求解 MHD 方程 组 存在 间断 的 黎 曼 问题 以 得 到 计算 区 域内 的 质量 、 能 量 、 动 量 等 各 
种 通 量 ， 从 而 利用 差分 格式 求解 . 图 9 给 出 , 4 CME 接近 1 AU 时 CME 和 扰动 太阳 风 的 
性 质 。 这 个 模型 用 来 探究 激 波 后 的 鞘 区 是 如 何 形成 的 ， 因 为 粒子 加 速 过 程 主要 依赖 于 这 个 
区 域 的 拓扑 结构 的 演化 情况 。 他 们 ”分 析 了 沿 着 磁力 线 的 密度 和 磁场 强度 的 跳跃 变化 ， 这 
些 量 反 映 了 等 离子 体 和 磁场 的 压缩 比 ,， 这 对 研究 扩散 加 速 、 一 阶 费 米 加 速 、 磁 泵 (magnetic 
pumping) 加 速 机 制 对 粒子 的 加 速效 果 具 有 重要 意义 。 他 们 还 认为 随 着 激 波 后 磁场 强度 的 突 
然 增 加 ,在 这 些 加 速 机 制 的 作用 下 , 粒子 的 能 量 能 够 达到 GeV HER. 
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3.2 AS-AkRRe 

A -H BRERA (CORHEL) 是 美国 集成 空间 天 气 建 模 中 心 (CISM) FAN. HAR 
和 日 球 层 两 个 模型 耦合 而 成 。 日 冕 模型 包括 三 种 选择 : (1) 用 势 场 源 表面 (PFSS) 模型 获得 
E ZRAK WSA (Wang-Sheeey-Arge) 模型 ”， (2) 能 量 方程 中 采用 多 方 过 程 的 MAS (the 
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MHD-Around-a-Sphere model) 模型 ; (3) 能 量 方程 中 采用 热 动 力学 过 程 的 MAS 模型 。 为 了 
计算 20 R, (Re: 太阳 半径 ) 或 30 Rs 至 1 AU 处 的 日 球 层 的 MHD 变量 , 日 球 层 模 型 可 以 选 
择 ENLIL 模型 或 者 日 球 MAS 模型 。 CORHEL 所 包含 的 模型 及 其 输入 输出 见 图 10 。WSA 
模型 以 光 球 视 向 磁场 观测 为 输入 , 用 PFSS 模型 获得 日 园 磁 场 , 根据 源 表 面 处 (通常 为 2.5 
R,) 的 日 蚁 磁场 相对 于 对 应 的 光 球 磁场 的 膨胀 因子 ,以 及 磁力 线 在 光 球 上 的 足 点 与 开 闭 场 边 
界 的 角 距离 ， 利 用 经 验 公式 确定 自 某 一 点 发 出 的 太阳 风 经 过 加 速 之 后 的 最 终 速度 。MAS 是 
一 个 三 维 时 变 MHD AIR A BRR, 其 中 磁感应 方程 用 Maxwell 方程 组 的 两 个 旋 度 方程 
和 欧姆 定律 代替 ， 引 入 矢量 势 表示 磁场 。ENLIL 模型 中 用 修正 的 TVDLF (Total Variation 
Diminishing Lax-Friedrich) 格式 计算 三 维 时 变 MHD 方程 。 此 模型 详 见 文献 [65,74] 。 
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图 10 ”CORHEL 模型 组 成 及 其 输入 输出 " 
红色 框 代表 输入 ， 蓝 色 框 为 相应 模型 , 绿色 框 代 表 输出 ， 未 来 潜在 的 发 展 内 容 用 “?” 表示 。 


3.3 ”太阳 行星 际 守恒 元 模型 

太阳 行星 际 守恒 元 解 元 (SIP-CESE) 模型 由 我 国 空间 天 气 学 国家 重点 实验 室 冯 学 尚 等 
Aw" 开发 , 是 从 太阳 表面 出 发 直至 行星 际 空间 的 三 维 数值 模型 。 控 制 方程 采用 守恒 形式 的 
MHD 方程 组 , 模式 的 特点 是 在 数值 格式 中 将 时 间 和 空间 不 加 区 别 作 为 一 个 整体 , 并且 对 空 
间 微 分 作为 待 解 变量 , 物理 量 在 解 元 (SE) 内 假定 是 光滑 可 微 的 , 而 在 相 邻 的 SE 之 间或 者 在 
守恒 元 (CE) 内 可 以 出 现 间断 。 该 方法 在 实际 计算 中 ，, 能够 有 效 避 免 日 地 空间 模拟 中 所 涉及 
的 球 壳 计 算 区 域 通常 产生 的 极 区 奇 性 和 网 格 收敛 性 问题 。SIP-CESE 模式 可 以 满足 磁场 散 度 
为 零 的 条 件 , 无 需 特别 处 理 , 并 引入 了 非 奇异 变换 ,使 得 自 适应 网 格 技术 (AMR) 容易 实现 ， 
计算 效率 得 以 提高 。 目前 该 模式 用 于 日 园 行 星际 太阳 风 背 景 的 研究 ”””、 太 阳 活 动 区 磁 
Rey? ”和 扰动 事件 的 模拟 7, 可 用 于 预报 和 研究 行星 际 背 景 参 数 和 日 时 物 质 抛 射 
及 其 驱动 激 波 传播 参数 。 冯 学 尚 等 人 ”用 改进 后 的 SIP-CESE MHD 模型 结合 AMR 技术， 
研究 了 4 个 不 同 卡 林 顿 自转 周 (CR) 内 (CR1967, CR2009, CR2060 和 CR2094) RRE 
景 ， 并 与 SOHO 和 OMNI 的 观测 结果 进行 对 比 , 图 11 是 计算 得 到 的 1 AU 准 稳 态 太阳 风 参 
数 与 OMNI 数据 的 对 比 , 所 得 到 的 太阳 风 密 度 、 速度 和 磁场 结构 呈现 出 该 时 期 背景 太阳 风 结 
构 的 典型 变化 特征 。 
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图 11 ”周期 内 的 太阳 风 参 数 在 1 AU 处 的 模拟 结果 与 观测 结果 的 对 比 
(a) CR 1967; (b) CR 2009; (c) CR 2060; (d) CR 2094. 红色 线 代表 模拟 的 结果 ; 绿色 线 代 表 OMNI 的 观测 数据 。 从 上 到 
下 依次 为 太阳 风速 度 人 、 等 离子 体 密度 n BET 和 径 向 磁场 B, 


3.4 “日 晃 行 星际 区 域 组 全 MHD (COIN-MHD) 模型 

EX RES X BA 6 (COIN) MHD 模型 ， 是 汉学 尚 等 人 对 开发 的 从 太阳 表面 出 发 直至 
行星 际 空间 的 又 一 个 三 维 数值 模型 ， 该 模式 将 从 日 面 到 行星 际 空间 的 太阳 风 计算 区 域 分 为 日 
园区 域 (大 约 为 1 ~ 22 R) 和 行星 际 空间 区 域 (18 Rs ~ 1 AU) 两 部 分 。 在 日 园区 域 , 先 用 
TVDLF 守恒 差分 格式 计算 出 新 时 间 步 的 值 , 再 将 磁场 部 分 用 MacCormack 了 [型 格式 的 计算 
结果 替换 掉 ; 而 在 行星 际 空间 区 域 用 MacCormack [格式 求解 与 时 间 无 关 的 定 态 MHD 方 
程 组 。 

ME AEA TE 2.5 R 处 构建 了 适合 不 同 太 阳 活 动 时 期 的 自治 源 表面 结构 , 为 原 
来 的 COIN-MHD 模式 提供 了 以 观测 为 基础 的 自治 初 边 值 条 件 。 此外， 他们 还 提出 一 种 新 的 
快捷 高 效 的 异步 并 行 时 间 推 进 法 , 并 应 用 到 三 维 太阳 风 及 CME 传播 的 MHD 数值 模拟 中 ; 
即 针对 日 晃 和 行星 际 空间 物理 性 质 的 不 同 , 根据 各 自 的 网 格 采用 不 同 的 计算 步 长 向 前 推进 ， 
得 到 了 精度 比较 高 的 结果 , 且 计 算 时 间 比 单一 步 长 缩短 了 8096 多 ,为 建立 快捷 的 太阳 风 模式 
提供 了 重要 参考 ” . 沈 芳 等 人 ”将 改进 后 的 模式 应 用 到 三 维 太 阳 风 背景 、CME 传播 及 相互 
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作用 等 方面 的 研究 中 , 以 2001 年 3 一 4 月 CME 相互 作用 事件 为 例 ， 比 较 了 两 个 CME 在 Ll 
点 (213 R,) 的 磁场 和 等 离子 体 参 数 的 模拟 结果 与 ACE 飞船 的 观测 结果 ( 见 图 12), 发 现 模 拟 
结果 可 以 再 现 和 解释 飞船 观测 到 的 多 重 磁 云 的 一 些 基 本 特征 。 

上 面 列 出 的 模型 中 , SWME 模型 以 强大 的 并 行 计算 能 力 为 基础 ， 可 以 实现 从 太阳 日 冕 到 
地 球 空间 环境 的 集成 研究 , 但 是 耗 时 长 。CORHEL 模型 中 对 日 园区 和 日 球 层 区 分 开 处 理 , 体 
现 了 模型 灵活 的 特点 。 其 中 MAS 可 以 实现 底部 边界 电磁 场 条 件 的 连续 更 新 ,因此 该 模型 可 
以 研究 日 冕 结构 的 演化 ,进而 为 日 球 层 模 型 提供 接近 实际 的 输入 ; 而 ENLIL 模型 中 采用 维 
度 分 裂 技术 提高 计算 效率 , 应 用 区 域 分 解法 实现 并 行 计算 。SIP-CESE FMEA BATE 
空间 采用 统一 的 处 理 方法 , 而 且 在 时 间 和 空间 也 采用 一 致 的 处 理 方法 , 它 区 别 于 所 有 现 有 太 
阳 风 模型 所 使 用 的 差分 格式 而 充分 体现 时 空 守恒 的 理念 。COIN 模型 针对 日 园 和 行星 际 的 不 
同 特点 , 采用 不 同 的 数值 格式 , 提高 了 计算 效率 ”。 上 述 模 型 所 得 到 的 结果 在 一 定 程度 上 能 
够 反映 观测 事实 , 在 模拟 太阳 风速 度 、 密度 、 总 磁场 的 变化 上 具有 较 好 的 准确 性 , 对 CME 或 
激 波 扰动 的 模拟 结果 能 够 定性 地 再 现 卫 星 实 地 观测 到 的 某 些 物理 量 的 变化 趋势 。 因 此 , 可 以 
研究 在 有 结构 太阳 风 缘 景 下 SEP 的 传播 ， 以 期 对 SEP 事件 有 更 加 深入 的 了 解 。 


4 BARRY SEP 事件 激 波 -粒子 模型 


Heras S A^ “首次 开发 了 关于 模拟 缓 变型 SEP 事件 的 激 波 -粒子 传播 模型 (Shock- 
and-Particle Model); Lario 等 人 " ”进一步 优化 了 Heras 的 模型 ， 并 提出 MHD 变量 与 粒子 
注入 率 的 Q(V R) 关系 。Rosa 在 Lario A" ^ 7. 建立 的 激 波 -粒子 模型 的 基础 上 进一步 完 
善 了 模拟 缓 变 型 SEP 事件 的 激 波 -粒子 模型 ， 这 个 模型 的 关键 在 于 加 入 了 观测 者 的 纬度 这 一 
因素 。 他 们 利用 Lario 等 人 "" 建立 的 MHD 变量 与 粒子 注入 率 的 Q(V R) 关系 , MZ MHD 
模拟 行星 际 激 波 , 首次 研究 了 激 波 强度 的 演化 、 激 波 加 速 粒子 注入 率 的 影响 、 质 子 通 量 剖 面 
的 变化 与 观测 者 所 在 纬度 的 可 能 相关 性 。Rosa A" "38 —48 MHD 模型 和 粒子 模型 结合 
在 一 起 , 提出 研究 缓 变型 SEP 事件 的 三 维 激 波 -粒子 混合 模型 , 其 显著 优势 是 可 以 通过 Lario 
等 人 提出 的 Q(VR) XAR, 将 粒子 模型 所 需要 的 注入 率 用 激 波 模 型 中 的 速度 计算 得 到 ， 从 
而 实现 了 对 粒子 模型 的 无 间断 输入 。 下 面 简 要 总 结 他 们 的 工作 , 本章 主 要 参考 Rosa 及 其 中 
相关 文献 的 内 容 。 

4.1 ” 激 波 模 型 

背景 太阳 风 初 始 条 件 采用 一 维 Parker"“ 流 体力 学 解 ， 磁 场 用 侦 极 子 场 ， 以 它们 作为 输 
A, 采用 时 间 松 弛 法 达到 稳定 状态 。 初始 的 磁 偶 极 子 位 形 消失 , 开放 和 闭合 磁力 线 区 域 形成 。 
为 了 形成 快慢 太阳 风 , 作者 采用 了 Groth A 的 半 经 验方 法 ， 即 在 能 量 方程 中 引入 额外 
的 加 热 项 。 轴 对 称 的 背景 太阳 风 模 型 用 三 维 球 坐标 (7,9,w) 下 的 MHD 方程 模拟 。 此 模型 的 
详细 介绍 参看 参考 文献 [96]。 

在 背景 太阳 风 中 从 加 一 个 高 密度 高 压 的 等 离子 体 团 来 形成 激 波 , 激 波 的 速度 vme 和 方 
向 Peme 给 定 。 扰 动 的 速度 和 密度 剖面 计算 公式 如 下 : 
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12 2001 年 L1 点 (213 Re) 处 磁场 和 等 离子 体 参数 的 模拟 结果 与 ACE 飞船 的 观测 结果 对 比 “” 

(a) 和 (b) 分 别 代表 ACE 的 测量 参数 和 模拟 结果 ,图 中 从 上 到 下 依次 为 : 地 心太 阳 黄 道 坐标 系 (GSE, geocentric solar 
ecliptic system) 下 的 磁场 强度 |B|(nT)，B, (nT)、 质子 数 密度 、 质子 温 度 、 速 度 以 及 等 离子 体 8。 其中, 蓝 色 实 线 表示 激 波 到 
达 L1 点 的 时 间 , 2 条 红色 虚线 代表 第 一 个 磁 云 的 开始 与 结束 时 间 , 绿色 虚线 表示 第 二 个 磁 云 的 开始 时 间 。 


a 一 Sone ( — ler) ; (6) 


其 中 , a 代表 密度 或 者 是 径 向 速度 ,aom。 是 a 的 最 大 值 , dom。 为 等 离子 体 团 的 半径 , d 为 球 
心 距离 。 

激 波 强度 用 VR= [Vr(d) —Vr(u)]/Vr(u) 表示 , u Al d 分 别 代表 激 波 的 上 游 和 下 游 。 此 
模型 中 , 粒子 的 加 速 过 程 被 忽略 , 激 波 被 处 理 成 一 个 源 项 。 
4.2 ”高 能 粒子 的 传输 模型 

在 研究 太阳 风 中 的 高 能 粒子 的 传输 时 需要 求解 聚焦 传输 方程 ，Rosa 等 人 采用 
Ruffolo ”发 展 的 聚焦 -扩散 方程 (见方 程 (7))， 其 中 t BN, p 为 投掷 角 的 余弦 值 ， 
是 日 心 距离 , p 是 粒子 的 动量 , v 是 粒子 的 速度 , y 是 径 向 方向 同 当地 磁场 方向 的 夹 角 , L 
是 磁场 聚焦 尺度 ，D,,, 是 扩散 系数 ，wsw 是 太阳 风 径 向 速度 。 此 方程 建立 在 混合 坐标 系 下 
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(r ts vw 是 固定 参考 系 下 的 变量 ; p, v, y 是 太阳 风 参 考 下 的 变量 ), 包括 粒子 沿 行星 际 磁场 
流动 , 投掷 角 扩 散 过 程 , 绝热 减速 ， 太阳风 对 流 和 磁场 聚焦 过 程 。 为 了 解释 组 变型 事件 ， 需 
要 在 运动 激 波 中 假定 一 个 连续 的 粒子 的 注入 , 在 方程 中 的 表现 形式 为 G(t, urp 它 与 相 空 
间 中 激 波 加 速 粒子 的 注入 率 Q 的 关系 为 G(t, pr, p) = A(7)Q(t, ur, p). Alr) 为 磁 通 量 管 的 
横 截面 , 在 cobpoint 位 置 上 计算 得 到 。Q 和 G 与 cobpoint 的 MHD 参数 相关 。 能 量 粒子 的 
传输 模型 假定 一 个 行星 际 磁场 为 Parker 螺 线 的 稳定 的 太阳 风 条 件 , 但 这 样 的 近似 只 有 当 上 
游 太 阳 风 没有 大 幅 扰动 的 情况 下 才 适 用 ; 而 且 假定 粒子 注入 发 生 在 cobpoint, 粒子 只 沿 着 磁 
力 线 方向 扩散 , 这 也 是 该 模型 的 不 足 之 处 。 对 于 粒子 的 传输 方程 ,Dalla 等 人 "” 用 单 粒子 一 
阶 绝热 理论 在 Parker 螺旋 场 中 验证 了 粒子 的 漂移 速度 在 SEP 传播 模拟 中 不 能 被 忽视 , 而 漂 
移 在 此 传输 方程 中 是 被 忽略 的 , 漂移 对 部 分 电离 的 SEP 重 离子 的 影响 相对 于 质子 来 说 更 强 。 
因此 , 粒子 的 传输 方程 也 需要 进一步 完善 。 


2 
DP ar) — cos yo x lv (: = ig Usw = Fe jsp) = 


TAE 
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43 ”主要 结果 

Rosa“” 分 别 在 0.4 AU 和 1.0 AU 处 假定 了 9 个 虚拟 的 观测 点 , 用 3D MHD 模拟 了 两 
个 行星 际 激 波 ( 快 激 波 和 慢 激 波 )， 运 用 激 波 -粒子 模型 ,合成 假定 事件 的 通 量 剖面 。 图 13 给 
出 了 快 激 波 条 件 下 由 粒子 传输 方程 模拟 得 出 的 质子 通 量 剖 面 ， 图 中 每 一 列 代 表 观 测 者 的 经 
E, 不 同 的 线 型 代表 观测 者 的 纬度 。 上 面 3 组 图 是 观测 者 处 于 0.4 AU 处 的 模拟 结果 , 下面 3 
组 图 是 观测 者 处 于 1.0 AU 的 模拟 结果 。 模拟 结果 显示 在 同时 考虑 径 向 距离 时 , 观测 者 纬度 
的 变化 可 能 导致 高 达 一 个 量 级 的 强度 变化 (比如 处 于 西 经 45° 的 观测 者 ), 说 明 在 缓 变型 事件 
的 模拟 中 ,纬度 与 峰值 强度 的 相关 性 是 不 容 忽 视 的 因素 ; 因此 , 缓 变型 SEP 事件 在 纬度 上 的 
演化 情况 和 经 度 方 向 上 的 演化 情况 一 样 也 要 得 到 研究 者 的 重视 。 


5 总 结 结 与 展望 


对 基于 MHD 太阳 风 模 拟 的 SEP 物理 模式 进行 研究 ， 关 键 是 建立 两 种 模型 的 联系 ,将 
描述 三 维 日 冕 物质 抛射 驱动 的 激 波 演化 模型 与 描述 粒子 在 行星 际 空间 中 加 速 和 传输 的 粒子 
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图 13 ”能 量 为 1 MeV( 蓝 色 线 ) 和 能 量 为 32 MeV( 红 色 线 ) 的 质子 通 量 剖面 ”” 
不 同 的 线 型 代表 观测 者 的 纬度 。(a) 观测 者 处 于 0.4 AU 处 的 模拟 结果 ; (b) 观测 者 处 于 1.0 AU 的 模拟 结果 。 


模型 耦合 起 来 。 在 之 前 模拟 激 波 的 模型 I^ ”中 , 一 般 假定 行星 际 磁场 是 Parker 场 , 假设 激 
波 面 为 平面 或 球面 ,而 且 不 考虑 粒子 在 下 游 的 传输 过 程 等 。 但 是 实际 情况 下 , 特别 是 在 太阳 
爆发 期 间 , 行星 际 的 磁场 结构 会 受到 剧烈 地 扰动 ,导致 粒子 在 行星 际 传播 时 ,加 速 和 传播 的 
过 程 复杂 很 多 。 因 此 ,未 来 的 研究 工作 中 可 以 尝试 加 入 真实 的 三 维 背景 太阳 风 模 型 和 三 维 
CME 激 波 模型 ， 并 且 考 虑 两 个 CME 在 很 短 的 时 间 尺 度 上 接连 爆发 时 对 激 波 加 速 粒子 的 影 
响 ， 从 而 更 加 深入 地 研究 粒子 的 加 速 和 传播 过 程 。 

对 于 组 变型 SEP 事件 的 混合 模拟 , 需要 结合 3 个 模块 : (1) CME 及 其 驱动 的 激 波 的 触 
发 和 传播 模型 ， (2) CME 驱动 的 激 波 的 粒子 加 速 模型 ，(3) SEP 在 三 维 行星 际 空间 中 的 传播 
模型 。 对 于 第 一 个 模块 , 建立 CME 驱动 的 三 维 激 波 从 太阳 附近 (比如 2 Rs) 到 行星 际 传播 
的 演化 模型 ， 以 此 获取 激 波 沿 磁力 线 的 cobpoint HBR, 作为 粒子 加 速 模型 的 输入 参数 。 对 
于 第 二 个 模块 , 建立 基于 聚焦 传输 方程 的 各 向 异性 的 一 维 激 波 加 速 模型 , 同时 考虑 沿 磁 力 线 


1 期 魏 稳 稳 ,等 : 基于 磁 流 体力 学 模拟 的 太阳 高 能 粒子 物理 模式 研究 进展 23 


的 激 波 加 速 和 存在 垂直 扩散 时 激 波 的 加 速 过 程 。 对 于 第 三 个 模块 , 可 以 尝试 加 入 实时 变化 的 
太阳 风速 度 及 接近 物理 真实 的 太阳 风 背 景 磁场 , 并 且 统 计 在 不 稳定 的 激 波 下 游 区 域 传 播 的 粒 
子 ,以 期 得 到 更 真实 的 粒子 在 三 维 行星 际 空间 的 传播 物理 过 程 。 综 上 所 述 ,建立 三 维 CME 
驱动 激 波 模 型 、 激 波 加 速 粒子 模型 和 SEP. 在 三 维 行星 际 空间 传播 模型 的 耦合 模型 ， 通 过 针 
对 具体 SEP. 事件 的 数值 模拟 , 将 观测 数据 分 析 工 作 与 模型 计算 密切 有 机 地 结合 , 不 断 检验 和 
改进 所 发 展 的 SEP 加 速 和 传播 模型 ， 是 未 来 研究 SEP 事件 的 具有 挑战 性 的 研究 方向 。 
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Abstract: Solar energetic particle (SEP) can destroy the spacecraft, and has great influence 
on human activities. Thus, it is of importance to forecast SEP events. SEP events can be 
divided by three general categories: impulsive SEP events, gradual SEP events and mixed 
SEP events. Impulsive SEP events, with characteristics of low intensity, short duration, and 
enhanced He-3, are correlated to solar flares. Gradual SEP events which become the main 
simulation branch of current research are related to coronal mass ejections (CMEs) and they 
usually have high flux in protons and last longer. Mixed SEP events maintain complex 
features of the former two kinds of events. 

Several theories and numerical models are applied to simulate gradual SEP events. Each 
model makes some assumptions to simplify the complex acceleration and transportation 
processes. No research to date has combined accurate shock evolution, particles injection 
and acceleration, with their interplanetary transport. In order to develop the physics-based 
SEP prediction model and improve the ability to forecast SEP events, it is necessary to self- 
consistently combine the magnetohydrodynamic (MHD) model that describes CME-driven 
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shock propagation with particle model that simulates the acceleration and transportation 
processes of particles. 

In this paper, firstly, we review current models, which have been used to describe 
gradual SEP events and their applications. Generally, there are two major approaches 
to model gradual SEP events: some studies include the acceleration mechanisms of SEPs 
induced by CME-driven shocks, while others assume a fixed particle injecting source at the 
shock. In addition, some researchers also consider the effects of perpendicular diffusion on 
SEPs propagation in three-dimensional interplanetary magnetic fields. These models can 
partially reproduce the observed properties for SEP events. Then, we make a brief review 
of numerical MHD simulation models, such as Space Weather Modeling Frame (SWMF), 
CORonal and HELiospheric (CORHEL) model, Solar InterPlanetary-Conservative Element 
Solution Element (SIP-CESE) model, COronal INterplanetary (COIN) model. All these 
models can be used to workout the propagation parameters of CME and CME-driven shocks, 
which are expected to provide inputs to the particle model. Finally, some discussions of the 


future work about how to combine MHD and particle models are presented. 


Key words: solar energetic particle; shock; MHD solar wind model; interplanetary transport 


